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Abstract: This is a study of holomorphic matrix models, the matrix models which 
underlie the conjecture of Dijkgraaf and Vafa. I first give a systematic description of the 
holomorphic one-matrix model. After discussing its convergence sectors, I show that 
certain puzzles related to its perturbative expansion admit a simple resolution in the 
holomorphic set-up. Constructing a 'complex' microcanonical ensemble, I check that 
the basic requirements of the conjecture (in particular, the special geometry relations 
involving chemical potentials) hold in the absence of the hermiticity constraint. I also 
show that planar solutions of the holomorphic model probe the entire moduli space 
of the associated algebraic curve. Finally, I give a brief discussion of holomorphic 
ADE models, focusing on the example of the A 2 quiver, for which I extract explicitly 
the relevant Riemann surface. In this case, use of the holomorphic model is crucial, 
since the Hermitian approach and its attending regularization would lead to a singular 
algebraic curve, thus contradicting the requirements of the conjecture. In particular, 
I show how an appropriate regularization of the holomorphic A 2 model produces the 
desired smooth Riemann surface in the limit when the regulator is removed, and that 
this limit can be described as a statistical ensemble of 'reduced' holomorphic models. 
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Introduction 

In the seminal paper [1], Dijkgraaf and Vafa proposed a beautiful conjecture relating 
matrix models to closed string theory on certain noncompact Calabi-Yau spaces. In 
its strongest form, this is meant as a relation between the partition function of a 
certain matrix model and the partition function of Kodaira-Spencer theory [21] on a 
dual noncompact Calabi-Yau threefold. The large N limit of this relation also leads to 
a matrix model description of corrections to the Veneziano-Yankielowicz potential of 
certain N — 1 supersymmetric field theories [1, 14] (moreover, the VY potential itself 
can be recovered from the non-perturbative part of the large N matrix integral). In fact, 
there are now independent derivations of the planar version of this conjecture by purely 
field-theoretic methods [2, 5, 4, 16], as well as first tests of it beyond the large N limit 
[6, 11]. More recent work on the subject can be found in [15, 8, 5, 17, 18, 13, 25, 26, 27]. 

As already pointed out in [1], a proper formulation of the conjecture should be 
given in terms of 'holomorphic matrix models', namely some version of matrix models 
involving 'contour' integrals in a space of complex matrices. This is very natural if 
one remembers that the conjecture was originally derived by considering the worldvol- 
ume theory of certain topological B-type branes, which is described by a reduction of 
holomorphic Chern-Simons theory [7]. Since the latter is formulated without reference 
to a metric, the proper description of the resulting matrix model should not involve 



2 



a hermiticity constraint. This expectation is borne out by the fact that the conjec- 
ture describes relations in the chiral ring of the dual field theory [4], and the latter is 
constrained by holomorphy. 

In [1], the authors chose the pragmatic approach of formulating the conjecture in 
terms of Hermitian matrix models, while pointing out that an appropriate holomorphic 
formulation should be given. The purpose of the present work is to initiate a systematic 
study of such holomorphic models, and show how they naturally tie up some loose ends 
encountered in various investigations and extensions of the conjecture. 

If naively taken at face value, the Hermitian formulation leads to various problems, 
some of which were already encountered (and overcome by a pragmatic prescription) 
in work aimed at testing it [6]. Among such issues, one can list the following: 

(a) Hermitian one-matrix models with generic polynomial potentials of odd degree 
(and complex coefficients) are ill-defined, since the real part of such potentials is not 
bounded from below along the real axis. However, such models are naturally required 
by the conjecture, which is supposed to apply without constraints on the degree of the 
potential. This issue was encountered in the work of [6], while performing a one-loop 
test of the conjecture for a cubic potential. As shown in Section 2 and Appendix 2, the 
pragmatic approach followed in [6] admits a natural justification in the holomorphic 
setup. 

(b) The large N spectral density ojq{z) = trj^ of a Hermitian one-matrix model 
(with matrix M) can only have cuts along the real axis. This means that the cuts of 
the hyperelliptic Riemann surface of [1] would be constrained to lie on the real line. In 
the Hermitian formulation, the conjecture would then imply that the dual Calabi-Yau 
space is constrained in a similar manner. Moreover, the Hermitian formulation leads 
to numerous problems in matching parameters and moduli spaces, since a polynomial 
of degree n always has n complex roots, but need not have n real roots. As clear from 
the work of [1] (and occasionally pointed out in subsequent work by the same authors), 
what is needed is a holomorphic extension of the Hermitian matrix model which would 
allow for an enlarged family of planar limits — namely, such a model should produce 
solutions whose large N eigenvalue distributions can be supported along arbitrary one- 
dimensional curve segments in the complex plane. 

In the present note, I show how these and related issues are solved by a direct 
analysis of holomorphic matrix models, whose construction follows a suggestion already 
made in [1]. In Sections 1 and 2, I explicitate the non-perturbative definition of such 
models 1 and study their convergence sectors, thereby refining one side of the conjecture. 

1 I should state from the outset that holomorphic matrix models are not the same as the so-called 
'complex matrix models' [19] sometimes encountered in the matrix model literature. Rather, they are 
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I also show (in Appendix 2) that issues like those encountered in [6] admit a natural 
resolution in this framework (which does recover the prescription used in that paper). 
Section 3 extracts the loop equations, equations of motion and the planar limit of such 
models, showing that most of the standard analysis employed for the Hermitian model 
carries through with certain modifications. Consideration of the large N limit leads to 
the algebraic curve of [1] , which is now unconstrained by any conditions on the location 
of the cuts. In fact, one can give a 'reconstruction theorem', which ensures that the 
holomorphic model probes the entire moduli space of this algebraic curve. This shows 
explicitly how such models solve the second issue listed above. Up to some details, 
the reconstruction result boils down to the well-known relation between the Riemann 
problem and singular integral equations. From this perspective, holomorphic matrix 
models give a sort of 'quantization' of the classical Riemann problem. 

The validity of the conjecture rests crucially on the special geometry relations men- 
tioned in [1], a more detailed justification of which was latter given in [9] (though only 
for the Hermitian case). To give a clear proof of such relations at the holomorphic 
level, Section 4 constructs a 'complex' microcanonical ensemble by introducing chem- 
ical potentials and performing a Legendre transform, thereby obtaining a generating 
functional which depends on averaged filling fractions; this allows one to show that 
the special geometry relations of [1] are a direct consequence of standard equations 
expressing chemical potentials in terms of the microcanonical generating function. In 
particular, these relations give the finite N analogue of the special geometry constraints. 

Section 5 considers holomorphic ADE models, focusing on a detailed analysis of 
the A 2 model. In this case, the problems one encounters by working naively with a 
Hermitian constraint are considerably more severe than in the one-matrix case. In 
fact, the Hermitian approach of [22] (with its attending regularization) would lead 
to a singular curve, which would always be constrained to have a certain number of 
double points. This would completely miss some of the gauge theory physics extracted 
in [23, 24] More precisely, such a constraint would force certain filling fractions to 
be identically zero, thereby contradicting the existence of the gaugino condensates 
obtained in [23]. As shown in Section 5, this issue is resolved quite naturally in the 
holomorphic A 2 model, by considering a regularization which is natural in that set-up. 
When removing the regulator, one obtains an ensemble of reduced holomorphic models, 
whose planar limit allows for non- vanishing values of all filling fractions. It is the large 
N limit of this ensemble which is relevant for the conjecture of [1]. 



a sort of 'square roots' of such models. 
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1. Construction of holomorphic one- matrix models and their 
eigenvalue representation 

Following a suggestion made in [1], we start by constructing holomorphic one- matrix 
models as multidimensional 'contour' integrals in a space of complex matrices. The 
construction is quite similar to that of a Hermitian model, except that the hermiticity 
constraint is replaced by a more general condition on the eigenvalues. After defining 
the model and studying its gauge-invariance, we extract an eigenvalue representation 
by choosing an appropriate multidimensional 'contour'. With this choice, one obtains 
an integral expression for the partition function which is formally identical to that of a 
Hermitian model, except that integration is performed over eigenvalues lying along an 
open path in the complex plane. 

1.1 The partition function 

Fixing a positive integer N, we let Mat^lC) denote the set of all N x N complex 
matrices. For any such matrix M, we let Pm(X) = det(XI — M) denote its characteristic 
polynomial. Define a subset M. of Mat n (C) as follows: 

M:= {M e Mat N (C) | p M (X) has distinct roots } . (1.1) 

This is an open submanifold of Maijv(C), consisting of matrices which are diagonaliz- 
able by general linear transformations: for every M in M, there exists an S E GL(N, C) 
such that: 

SMS" 1 — D :— diag(A! . . . \ N ) , (1.2) 

where Xj are the roots of Pm(X). The later coincide with the eigenvalues of M, so that 
its spectrum is given by: 

a(M) = {X 1 ...X N } . (1.3) 

Consider a connected, noncompact and boundary-less submanifold r of M. of real 
dimension equal to the complex dimension of M.: 

dim n T = dim c M = dim c Mat N (C) = N 2 . (1.4) 

Also consider the standard holomorphic symplectic form on Mat N (C): 

w = AijdMij . (1.5) 

The sign of the right hand side of course depends on the ordering of pairs € 
{1 . . . N} x {1 . . . iV}, and we shall use the lexicographic ordering: 

w = dM u A ... A dM 1N A dM 2l A ... A dM 2N A ... A dM Nl A ... A dM NN . (1.6) 
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This is implicit in all such expressions encountered below. 

Fixing a polynomial W(z) = X^ =0 £ m z m with complex coefficients, we define the 
holomorphic one-matrix model by the partition function: 

Z N (T,t) = ± J r we-»» W M , (1.7) 
where Af is a normalization factor to be fixed below. 

1.2 Gauge invariance 

Consider the GL(N, C) action r on Mat jv(C) given by similarity transformations: 

r(5)M := , 5 G GL(N, C) . (1.8) 

This clearly stabilizes At; in particular, note that the characteristic polynomial of M 
is r-invariant: 

Pt(s)m(A) =Pm(A) . (1.9) 

The action S(M) = NtrW(M) of our model is obviously invariant with respect 
to (1.8). The same is true of the holomorphic measure w: since [r(S)M]ij = a^^M^, 
with aij t ki = SikiS -1 )^ (i.e. a = S (8) (5' _1 ) T ), we have det(a) = 1 and: 

t(S)*(w) =w . (1.10) 

It follows that the model admits the GL(N,C) gauge-invariance (1.8), provided that 
one chooses the integration manifold T such that it is stabilized by the gauge-group 
action. 

Observation Choosing S to be a permutation matrix (i.e. Sy = with a 

a permutation on N elements), equation (1.10) shows that the holomorphic measure 
(1.5) is invariant under joint permutations of indices: 

dM n A ... A dM 1N A dM 21 A ... A dM 2N A ... A dM N1 A ... A dM NN = (1.11) 
dM a{1)a(1) A ... A dM a(1)<j{N) A dM a(2 )a(i) A ... A dM^^aiN) A ... A dM <r ( JV ) (7 (i) A ... A dM a(N)cr{N) , 

a relation which can also be checked directly. 

1.3 Eigenvalue representation 

One can derive an eigenvalue representation of (1.7) as follows. Let 7 : R — > C be 
an open curve in the complex plane, which we assume to be immersed and without 
self- intersections. We shall take Y to be the following subset of M.: 

r( 7 ) := {M G M I a(M) C 7} , (1-12) 
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which obviously satisfies (1.4). We denote Z^{Y{^),t) by 2^(7, i). Then the integral 
over gauge-group orbits can be performed as explained in Appendix 1, with the result: 

Z N (l,t) = ^(-l) N2 ^ 2 ^ hvol(H)Z N ( 1 ^ , (1.13) 

where hvol(H) is the 'holomorphic volume' of the complex homogeneous space H = 
GL(N, C)/(C*) N (see Appendix 1) and 

Z N (>y,t) = [d\ 1 . . . [ dX N H (\ - Xj)e- N ^U w ^ . (1.14) 

is the eigenvalue representation of our model. The holomorphic volume hvol(H) will be 
discarded (together with the sign prefactors) by choosing A/" = (— l)-^ 2 ^ -1 )/ 2 -^ hvol(H) 
in (1.7). Expression (1.14) is formally identical with the eigenvalue representation of 
the Hermitian matrix model, except that the eigenvalue integral is performed along the 
contour 7 in the complex plane. The pragmatic reader can take (1.14) as the definition 
of our model. 



2. Convergence sectors 

It is clear from expression (1.14) that convergence of our partition function depends on 
the choice of 7. In this section, we shall characterize the 'good' choices of 7 in terms 
of certain asymptotic sectors of the model, described in terms of cones partitioning 
the complex plane. Such cones can be identified by performing an elementary analysis 
of the behavior of the integrand. As we shall see below, this allows us to make non- 
perturbative sense of models with polynomial potentials of odd degree. 

To extract the relevant behavior, let z = re ld with r > and 9 G R/27rZ and 
tj = rje j with r,- > and 9j G R/27rZ. Then the potential takes the form: 

n 00 

w(z) = t + Y, r f j cos U - 6 j) + i 2 r f j sin C^ - Oj) . (2- 1 ) 

and the behavior of |e -ArH/(z )| = e ~ NReW ( z ) for r — > 00 is dictated by the contribution 
rj r JO cos (jo — 9j ), where j = jo (9) is the largest j such that rjCOs(j9 — 9j) 7^ 0. 
Namely, e~ NW ( z ) is exponentially decreasing/increasing depending on whether cos(jo# — 
9j ) is positive /negative. In the very non-generic case when rj cos(j9 — 9j) vanishes for 
all j — 1 . . . n, the quantity e~ NW ^ oscillates indefinitely as r — > 00. 

Let us fix 9 such that cos(n9 — 9 n ) 7^ 0. Then e ~ NW ( z ) [ s exponentially decreasing 
as r — > 00 if and only if cos(n9 — 9 n ) > 0, which gives: 

9= a + 0n + 7T— with fc = 0...n-l and a G (-vr/2, tt/2) . (2.2) 
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This defines n angular sectors (=open cones with apex at the origin) in the complex 
plane, which we denote by A k , k — 0. . .n — 1. We also define complementary sectors 
B k through: 

6 = Q + 6n +Tr 2k+ 1 with fc = 0...n-l and a G (-tt/2, tt/2) ; (2.3) 
n n 

these are the sectors for which cos(n6> — n ) < 0. The sectors A fc and B k arise consecu- 
tively with respect to the trigonometric order and are separated by rays originating at 

2 = 0. 

To make the integral (1.14) absolutely convergent, we shall require that the curve 
7 asymptotes to some straight lines d±(t) = ±v±t + fi± for t — > ±oo, such that the 
corresponding asymptotes lie inside two of the convergence sectors A k . That is, we 
require: 

3 Una ("f(t) =F v±t — n±) = (2.4) 

t^±oo 

3 lim i(t) := ±v± . 

for some u± E A k± and some /a± £ C. Here ^'(t) = and k± E {0 . . . n — 1}. Since 
the integrand of (1.14) is holomorphic, it immediately follows that the partition function 
is independent of the choice of 7 as long as this contour has asymptotic behavior (2.4) 
with v± belonging to fixed sectors A k± . In particular, the integral does not depend 
on the precise values of v± and /i±, but only on the convergence sectors A k+ and A k _. 
Therefore, we obtain n 2 possible values Z(k-,k + ,t), indexed by the sectors A k _ and 
A k+ (figure 1). A complete definition of the model requires that we specify one of the 
'phases' k + ), together with the potential W. 

It is also clear that Z{k_ 1 k + , t) vanishes if k + = k_, and that we have the relation: 

Z(k_,k + ,t) = (-l) N Z(k + ,k_,t) , (2.5) 

which results upon reversing the orientation of 7. Therefore, it suffices to consider the 
n(n — l)/2 'phases' indexed by pairs k+) with k- > k + . 

2.1 Scaling symmetry 

Let q be a non- vanishing complex number. Writing z :— |, we have: 

W(z) = W q (x) , (2.6) 

where W q (x) = J2]=o t^x j , with '■= ^j- Performing the change of variables A := | 
in (1.14) gives: 

Z N {n,t) = ± s Z N ( lqi tto) , (2.7) 



8 




Figure 1: Convergence sectors of the holomorphic matrix model. We show the case n = 
degW = 4, with 64 = and a contour belonging to the sector (£;_, k + ) = (1, 0). 

where 7 g is the path defined through: 

7 ,(t) = <n(t) (2.8) 

for all real t. 

For the choice q n := rl/ n e~~^ , the change of variable 2 := — gives £^™) = 1, so 
the transformed potential W qn has unit leading coefficient. Hence the model depends 
'trivially' on the parameter t n , and we can set t n = 1 and 8 n = without loss of 
generality. 

Also note that choosing q := a k with a = and k an integer allows one to rotate 
7 by any multiple of — . Since this does not change the convergence sectors (because 
t^ k ^ = t n ), it leads to the relation: 

Z N (k + k_,k + k + ,{t m }) = a N ' 2k Z N (k„,k + ,{a mk t m }) . (2.9) 

Hence it suffices to consider the 'phases' (fc_, 0) with k- — 1 . . .n — 1. 

Observation: If one increases the degree of W, then the convergence sectors be- 
come progressively narrower. Allowing for the case n = 00 (i.e. for an entire func- 
tion W) is often a useful device in the theory of matrix models (the best known 
example is the matrix model of Kontsevich [12] and its generalizations). In this 
case, the convergence structure of the holomorphic model can be quite different from 
that discussed above, and depends on the precise asymptotic behavior of the entire 
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function W. A simple example is provided by the choice W{z) = e z , which gives 
| e -AW(z)| — e -NRe w(z) _ e -Ne x cos(y) ^ w ;h ere z — x _|_ iy w jth x,y real. Then the 

convergence sectors are horizontal strips defined by the condition: 

cos(y) > y e (~ + 27rk, | + 2vrfc) , k G Z . (2.10) 

Correspondingly, we obtain a convergent partition function by taking 7 to satisfy (2.4), 
where now v± > and /i± belong to two such bands (figure 2). This example should 
serve as warning against the idea that one can recover models with power series poten- 
tials as naive limits of polynomial models. 



571-/2 






7T/2 


d < 




-7T/2 




-37T/2 

-5tt/2 





Figure 2: Convergence sectors for an entire potential W = e z , and a good choice of contour 
for such a model. The filled-in regions are forbidden sectors for 



2.2 Even and odd degree potentials, and the Hermitian matrix model 
2.2.1 The case n =even 

In the case n =even, the convergence sectors of the model appear in pairs A k , A k+n / 2 
which are symmetric with respect to the inversion z —>■ —z. By the discussion above, 
we can take t n positive without loss of generality. Then 9 n = and the two sectors 
A and A n / 2 lies opposite to each other and are bisected by the positive real axis. 
Picking k_ = n/2 and k + = 0, we can describe the 'phase' Z(n/2,0,t) by taking the 
curve 7 to coincide with the real line (with its usual orientation). Then the partition 
function (1.14) reduces to the usual eigenvalue representation of the Hermitian one- 
matrix model. This gives a partial justification for the formulation used in [1]. 
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2.2.2 The case n =odd 

In this case, one has an odd number of convergence sectors Ak, whose images under the 
reflection z — > —z are the 'bad' sectors B k+ [ n / 2 ]- Picking t n to be positive as above, it is 
clear that the matrix model cannot be defined by choosing 7 = R, since the associated 
integral would diverge; this amounts to the basic observation that the Hermitian matrix 
model is ill-defined for generic complex polynomial potentials of odd degree. However, 
one can certainly make sense of the holomorphic model by working with any of the 
good 'phases' (k-, k + ) — one simply picks a contour 7 such that u + and z/_ belong to 
some of the convergence sectors A^. 

The fact that the Hermitian model cannot be relevant for odd degree potentials 
is related to problems observed in [6], which in that paper were avoided by declaring 
that certain matrices occurring from a perturbative analysis should be anti- Hermitian. 
A systematic formulation of this idea is to consider the holomorphic matrix model 
instead. To substantiate this proposal, Appendix 2 shows that the procedure of [6] 
admits a simple justification in the framework of holomorphic models. 

3. Loop equations, equations of motion and the large N limit 

In this section, we study the loop equations and equations of motion of the holomorphic 
model, as well as its planar limit. As we shall see below, much of the standard fare 
of Hermitian models can be extended quite directly to the holomorphic case (though 
there are a few modifications). The most important novel feature of holomorphic matrix 
models is that they can explore an enlarged set of planar limits, thus probing the entire 
moduli space of a certain family of algebraic curves. This fact, which is essential 
for the conjecture of [1], will be established explicitly by proving a reconstruction 
theorem which associates a planar solution of the model with an arbitrary algebraic 
curve belonging to this family. 

Let ds be the length element on 7 and s the length coordinate along this curve, 
centered at some point on 7. We let A(s) denote the parameterization of 7 with respect 
to this coordinate and write Aj = A(sj) accordingly. Mimicking usual constructions, we 
define the spectral density by: 

1 N 

, (3-1) 

JV 3=1 

where 5(s — Sj) is the delta- function in the coordinate s (equivalently, the 5-function 
along 7 with respect to the measure induced by the length element). Note the normal- 
ization condition: 

/oo 
dsp{s) = 1 . (3.2) 
-00 
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Also consider the resolvent of M: 



R(z) : = 



1 



and its normalized trace: 



z-M ' 



1 x 1 1 1^1 f 

u(z) = —trR(z) = —tr — = — V = / 

v; iV v; N z-M N^z-Xi J 

1=1 1 



(3.3) 



ds- 



P{s) 



z-X(s) 

In the following, we shall need the Sokhotsky-Plemelj formulae: 

lim / dX'- \ = V f dX'—?— Ti* for A G 7 , 

e^o+ U A — A' ± en A U A — A' 



(3.4) 



± era(A) 

which we also write symbolically as: 

1 



(3.5) 



lim 



V- 



1 in S( /\ 



e ^o+ X{s) - X(s') ± en(s) X(s) - X(s') X(s) 



(3.6) 



Here V stands for the principal value and n(s) = iX(s) is the normal vector field to 7 
(figure 3). 
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Figure 3: The normal vector field n(s) = it(s), and the tangent vector field t(s) = A(s) of 
7. Note that \n(s)\ = \t(s)\ = 1 since s is the length coordinate along 7. 



3.1 Loop equations 

To extract the loop equations, we follow the method of [10]. For this, start with the 
identity: 



[dx 1 ...[dx N {:^- n (a* - AOe-^- w ^ ] — 

J ~f J ~' i=l OA i k+l Z ~ 



A, 



. 



(3.7) 
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Performing the partial differentiation, we write this as: 

K 1 j\rJ2 W '( X i) | 1 ) = o . (3.8) 

Using the decomposition: 

f 1 



(z — a)(z — (3) a — (3 



z — a z — f3 



(3.9) 



to simplify the last term and combining the result with the first gives: 

Defining the polynomial: 

1 f W\z) - W'(Xj) r W'(z)-W'(X(s)) 

relation (3.10) gives the following form of the loop equations: 

(u(z) 2 )-W'(z)(u(z)) + {f(z))=0 . (3.12) 
3.2 Equations of motion 

Writing A; = A(sj), the partition function (1.14) becomes: 

Z N ( r t) = Jd Sl ...J d ajv n^( a *)n( A ( a *)-A(a i ))e- Jv ^=i wr W'')) 

1 = 1 iy^j 

= j 'ds 1 ... J ds N e- NS ^f {si - SN) , (3.13) 
where A(s) := and: 

N i i TV 

S eff ( Sl ...s N ) = J2 W(X( Si )) - - £ H\ - A,-) - T7 E ln ^ • (3.14) 

j=l JV i# iv 1=1 

Extremizing with respect to gives the equations of motion: 

The prime means that the sum is taken only over j ^ i. The last term is a curvature- 
induced contribution which is sub leading in 1/JV. It is also easy to check that the 
equations of motion imply: 

a**)' - W^Mz) + /(,) + + _ g ^ = . (3.16) 

This 'operator' equation holds only for solutions of (3.15), unlike the Ward identity 
(3.12). 
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3.3 The large N limit 

For any quantity <j), consider the large N expansion: 



oo 



<*> = £ 



(3.17) 



3=0 



with coefficients 4> , (fix etc. In particular, we have (p(s)) = p (s)+O(l/N) and (uj(z)) = 
ujq(z) + 0(1/N). In the large N limit, the eigenvalues Xj are replaced by the planar 
spectral density po(X) supported on the curve 7. Note that this quantity is complex, 
since its definition involves averaging with respect to the complex integrand of (1.14). 
As usual, we have 



so that uj becomes an analytic function whose cuts I a are forced to lie on the curve 7 
(figure 4). We shall assume for simplicity that all cuts are of finite length. 



Figure 4: Cuts of ljq. We also show a closed contour T surrounding all cuts and a point z 
in its exterior. 

In the planar limit, the average of a product can be replaced by the product of 
averages and the loop equations (3.12) reduce to the algebraic constraint: 




(3.18) 




z 



r 




W'(z)u o (z)+f o (z)=0 , 



(3.19) 



where 




(3.20) 
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is a polynomial of degree n — 2 with complex coefficients. Since p(s) is normalized by 
(3.2), equation (3.20) shows that the leading coefficient of fo(z) equals nt n . 
In the planar limit, the equations of motion (3.15) become: 

The curvature term in (3.15) drops out, since it is subleading in 1/N. The algebraic 
constraint (3.19) can also be obtained from (3.21) by standard manipulations. Writing: 

u (z)=uo(z) + ^W\z) , (3.22) 
relation (3.19) shows that u (z) is one of the branches of the planar affine curve: 

u 2 -\w'(z) 2 + Mz) = . (3.23) 

Since the branch cuts of ujq (and thus of uo) must lie along 7, it is clear that the 
polynomial /o in equation (3.23) is constrained by the choice of this curve. However, 
changing 7 without changing its asymptotes allows one to describe any position of the 
cuts in the complex plane, as long as all these cuts have finite length. This effectively 
eliminates the constraints that would be present in the Hermitian case (for which 7 
would be forced to coincide with the real axis). 

3.4 Reconstruction of a planar solution from the Riemann surface 

Given the algebraic curve (3.23), we now show how one can use it to recover an appro- 
priate 7 supporting a spectral density po satisfying (3.2), (3.18), (3.20) and the planar 
equations of motion (3.21). This proves that the holomorphic matrix model is free to 
explore the whole relevant piece of the moduli space of (3.23), unlike the Hermitian 
matrix model. In the planar limit, the holomorphic model reduces to the singular in- 
tegral equation (3.21). Up to some minor details, the existence of a one to one relation 
between solutions of this equation and members of a family of Riemann surfaces boils 
down to the well-known relation between the Riemann problem and singular integral 
equations. 

To see this explicitly, assume that one is given a complex degree n — 2 polynomial 
fo(z), subject only to the constraint that its leading coefficient equals nt n . Consider 
the associated curve (3.23). Denoting its branch points by a a , b a (with a — 1 . . . n — 1), 
choose branch-cuts I a connecting a a and b a (note that we allow I a to be curved). This 
defines two determinations, which we call uq and u\ = —u$. More precisely, uq is the 
determination which behaves as ~W'(z) for large z. Let us define u by relation (3.22) 
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and choose a curve 7 such that I a C 7 for all a. We let s be its length coordinate and 
A = A(s) the associated parameterization. Choosing the normal n(s) = iX(s) (figure 3) 
we define: 

Po(s) := X(s) lim [u (\(s) - en(s)) - u (X(s) + en(s))} (3.24) 

for every s in 7. Then p vanishes outside of I a and (3.24) and the Sokhotsky formulae 
(3.5) imply: 

/*-^=/£— • (3 - 25) 

7 2: — A(s) Jr 2m 2 — x 

where T is a contour surrounding all cuts but not the point z (see figure 4) and the 
last equality follows by deforming this contour toward infinity to pick the contribution 
from the residue at x — z. This shows that relation (3.18) holds. Using (3.23) and 
the fact that the leading coefficient of /o equals nt n shows that tv(z) = 1/z + 0(1/ z 2 ) 
for large \z\, and combining this with (3.18) shows that p satisfies the normalization 
condition (3.2). 

Since the cuts I a connect the branches uq and —u , we have: 

lim uq(\ + en) = lim [— uniX — en)] , A e I a (3.26) 

so that lim e ^ + [cj (A + en) + uj (X — en)} = W'(X) for A G I a . Combining this with 
equation (3.18) and using the Sokhotsky formulae (3.5) shows that po satisfies the 
planar equations of motion (3.21) along the cuts. 

To prove (3.20), we use relation (3.18) to compute: 



(3.27) 

where we used the identity (3.9) and symmetry of the resulting integrand with respect 
to the substitution s <-> s'. Using (3.21) in the right hand side of (3.27) gives: 

= J^«M^§ = - j dsK(.) W 'M I |p + W(z»Mz) . (3.28) 

Comparing with (3.23) and (3.22) shows that equation (3.20) holds. 

This construction produces a planar solution of the holomorphic model which re- 
covers any curve of the form (3.23), for an arbitrary choice of degree n — 2 polynomial 
fo(z) with leading coefficient nt n . Unlike the Hermitian model, the holomorphic model 
explores the entire family of curves (3.23) in its planar limit. This is a basic pre-requisite 
for the conjecture of [1]. 
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It is also clear that the precise choice of cuts I a is irrelevant, as long as they 
connect the given pairs of branch points a a ,b a . Thus one can use any 2 curve 7 in this 
construction, provided that it passes through all the branch points of (3.23). 

4. The microcanonical ensemble 

The framework of [1] requires that the model obeys certain filling fraction constraints. 
In this section, I explain how one can impose such constraints on the finite N model 3 . 
The relevant conditions are easiest to formulate by employing a microcanonical en- 
semble. As we shall see below, the original path integral defines a (grand-) canonical 
ensemble at zero chemical potentials. This allows one to recover the microcanonical 
generating function by introducing chemical potentials (which are canonically conju- 
gate to the filling fractions) and then performing a Legendre transform to replace the 
former by the latter. 

4.1 The (grand-) canonical partition function associated with a collection of 
domains 

We start by fixing a collection of domains D a (a = 1 ... r) in the complex plane, chosen 
such that their interiors are mutually disjoint and such that: 

uUA* = C . (4.1) 

We shall assume that 7 intersects each closure D a along a single curve segment A a , 
where A a lie in ascending order on 7 with respect to its orientation (figure 5). Condition 
(4.1) implies: 

U' Q= A* = 7 . (4.2) 

For simplicity, we shall take D a to be infinite strips arranged as in figure 5 . 

Letting Xa denote the characteristic function of D a , consider the matrices Xa(M) 
defined by holomorphic functional calculus: 

2 Nondcgenerate and without self-intersections. 

A clear formulation of such constraints beyond the strict large ./V limit is important in studies 
of SO(N) and Sp(N) matrix models, whose large N expansion contains terms of order 1/N which 
must be interpreted as contributions to the dual field theory (rather than gravitational contributions, 
which start at order 1/N 2 ). The formulation given below allows one to give clear proofs of relations 
between the strictly planar and 0(1 /N) contributions to the (microcanonical) partition function, 
thereby strengthening arguments like those given in [8]. Similar problems are encountered in more 
general orientifold models [28]. 
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Figure 5: A choice of strip domains in the complex plane. We also show two of the bounding 
contours (namely 1^ and T§). In this example, we take the domains to be infinite strips; this 
assures us that 7 and its deformations will cut each such domain along a non-void interval. 

with T a a (counterclockwise) contour bounding D a in the complex plane (figure 5). 
The matrix Xa(M) equals the projector on the space spanned by those eigenvectors of 
M whose eigenvalues lie in D a . 

Let f a := jftrXa(M) = — Y,f=i Xa(Xj) be the filling fractions for the domains D a . 
One has: 

f a = [ dsp(s) Xa (\(s)) = I ^co(z) , (4.4) 

since u(z) has simple poles with residue j? at each eigenvalue Xj. Relation (4.2) implies 
Eq=i Xa(^) — 1) so the filling fractions are subject to the constraint: 

E fa = 1 • (4-5) 

Q = l 

Picking complex chemical potentials fi a , we consider the (grand-) canonical ensemble 
associated with our collection of domains: 

z N ( r ,t,v) = 1 / dMeM-^w-^EL.^^w] 

N Jv(i) 

= I dX 1 ...[ ^n^-A^e^^-^^ 2 ^^ , (4.6) 
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which is an analytic function of \i a . The original partition function results by setting 
\i a = 0, and thus it corresponds to a (grand-) canonical ensemble at zero chemical 
potentials. Notice that the (grand-) canonical partition function can be written: 

z N (r,t^)= £ N ^\ [ e- N ^^z Nl ... Nr ( T ,t) , (4.7) 

N 1 +...+N r =N iVl iVr ' 

where: 

Z Nl ...N r (r, t) = [ d\i . . . ( d\ Nl . . . I dX Nl+ ... +Nr _ 1+1 ... I dX N J| (A* - Xj)e~ N ^^ W{ - Xj) 
j 1 j 7 j 7 j 7 

(4.8) 

Observation One can consider a more general version of microcanonical ensemble 
based on a finite partition of unity, i.e. a finite collection of smooth complex-valued 
functions <f> a (z,z) satisfying the constraint: 

E^ = 1 • ( 4 - 9 ) 

a 

In certain ways, this is preferable to the approach taken above, since it may be techni- 
cally desirable to avoid having to deal with characteristic functions. The entire discus- 
sion of this section extends easily to this more general set-up. 
Introducing the (grand-) canonical generating function: 

F N (r,t,n) := -^\nZ N (r,t,/i) , (4.10) 
we have the standard relation: 

£? = M ■ < 4 - n > 

Here and below, the brackets (. . .) denote the expectation value taken in the (grand- 
)canonical ensemble. 

4.2 The microcanonical generating function 

Following standard statistical mechanics procedure, we define: 

S a := tt~ J~ (4.12) 
and perform a Legendre transform to extract the microcanonical generating function: 
F( 7 , t, S) := 5 a // Q ( 7 , t, S) - ^( 7 , t, n a (t, S)) , (4.13) 
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which is an analytic function of the complex variables t a and S a . In this relation, [i a 
are expressed in terms of t and S by solving equations (4.12). The constraint (4.5) and 
equation (4.11) show that S a satisfy: 

N 

£s Q = i » ( 4 - 14 ) 

a=l 

so we can take Si . . . £y-i to be the independent variables. Then equations (4.12) 
express /i a as functions of tj and these coordinates, and equation (4.13) implies: 

dF 

Ha ~ = ^tt for a = l...r- 1 . (4.15) 

This gives the chemical potentials as functions of t and S a . Note that fi a are only 
determined up to a common constant shift; this is due to the constraint (4.14) on S a . 
Working with F(^y, t, S) amounts to fixing the expectation values of the filling fractions 
by imposing the quantum constraint (4.12): 

{fa) = L^n {uj{z)) = Sa ' (4 ' 16) 
with S a treated as fixed parameters. We stress that this condition is only imposed on 
the expectation values of the filling fractions. 

4.3 Chemical potentials at large N 

We will now show that the large N chemical potentials can be expressed as B-type 
periods of the algebraic curve (3.23), thereby proving that the special geometry relations 
of [1] hold at the holomorphic matrix level. The argument below is an adaptation of 
that given in [9] 4 , combined with the definition of the microcanonical ensemble given 
above. In particular, we show that the special geometry relations are simply the large 
N limit of the standard equations (4.15). Hence the chemical potentials fJ>o> £1X6 the 
'quantum' (i.e. finite N) analogues of the B-type periods of [1], while the averaged 
filling fractions are the 'quantum' analogues of the A-type periods. This captures the 
beautiful intuition of [1]. 

For this, we start from the expression of the (grand-) canonical generating function 
in the planar limit: 

Fo(l,t,n) = I dsW{\{s))p Q {s)-V fds f ds'K(\(s),\(s'))p (s)p (s') + j2fi a I ds Po (s) 

Q=l ^ a 

(4-17) 

4 In [1], Dijkgraaf and Vafa gave a beautiful intuitive justification for the special geometry relation 
n a = §§ s - + const, where F a is the planar limit of the (microcanonical) generating function while S a 
and I1 Q are identified with periods of the curve (3.23). A derivation of this relation (in the context of 
the Hermitian model) was later given in [9], upon using older results of [20]. 
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where: 

K(\,\') :=ln(A-A') . (4.18) 

In what follows, we shall assume that each cut I a lies inside a corresponding domain 
D a ; in particular, we assume that the number of cuts coincides with the number of 
domains. With this assumption, we have S a = (f a ) = Jj dspo(s) and the Legendre 
transform of (4.17) gives the planar limit of the microcanonical generating function: 

F ( 7 , t,S)=vJdsJ ds'K(X(s), X(s'))p ( S )p (s') - J dsW(X( S ))p (s) , (4.19) 

with the constraints: 

[ dsp {s) = S a for a = l...r . (4.20) 

Remember that we allow I a to be curved intervals connecting the branch points a a and 
b a of the algebraic curve (3.23). In the generic case, none of the cuts is reduced to a 
double point and one can take r — n — 1. 

4.3.1 The primitive of u along 7 

To extract the large iV chemical potentials, consider the 'restriction' of u along 7, 
which we define by: 

wg(s) := l lim [u (X(s) + en(s)) + u (X(s) - en(s))} . (4.21) 
Here n(s) = iX(s). 

If A(s) is a point of 7 lying outside the union of I a , then Uq(s) equals Uq(X(s)), the 
quantity obtained by substituting X(s) for z in the expression: 

u (z) = cvo(z) - \w\z) = J ds'^^ - \w'{z) . (4.22) 

Using (4.22) in the definition of u p gives: 

ul{s) :=VJ ds'p (s')K(X(s), X(s')) - ±W'(\( 8 )) , (4.23) 



where: 

K ^^ = Y^x (424) 

is the integral kernel appearing in the planar loop equations and where we used the 
Sokhotsky identities (3.5). 

Consider now the function : 7 — > C defined though: 

0(s) := 2V J ds'K{X( S ),X( S '))p (s') - W(X(s)) . (4.25) 
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Noticing that K(X, A') = ^K(X, A'), we have: 

±<P(s) = 2\(sH(s) . (4.26) 

As clear from (4.23), the planar equations of motion (3.21) amount to the require- 
ment that Uq vanishes along each of the curve segments I a : 

u p (s) = for X(s) G / := U r a=1 I a . (4.27) 

This means that is constant along each of these intervals: 

4>{s) = i a = constant on I a , (4.28) 

The jump in the value of u between consecutive cuts can be obtained by integrating 
(4.26): 

- 6* = 2 dXuo(X) , (4.29) 

Jb a 

where we used d\ = X(s)ds. This integral is of course taken along 7. 
4.3.2 The planar chemical potentials 

Differentiating (4.19) with respect to S a for some a < r and using relation (4.25) gives: 

^ 0) - ^ = 7^o(7, t,S)= f ds^-^s) = - , (4.30) 

where fj,^ are the planar limits of the chemical potentials. To arrive at the last equality, 
we used equation (4.28) and noticed that J I/3 8 ^f A ) = tJ- J T dsp (s) equals 5 af3 if (3 < r 
and —lUP — r (by virtue of (4.20) and (4.14)). Relation (4.30) shows that the planar 
chemical potentials coincide with the quantities £ a , up to a common additive constant. 
Using relation (4.29), we obtain: 

- ^ = 2 dXu (X) (4.31) 

Jb a 

As explained above, the quantity u (z) has cuts precisely along the intervals I a . 
Since the other branch of (3.23) is given by U\{z) = —u (z), this allows us to write 
(4.31) in the form: 

^ - = dz[ Ul (z) - u (z)\ = I dzu(z) , (4.32) 
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Figure 6: Choice of A and B-cycles on the large N Riemann surface (for convenience, 
we represent the curve 7 as a straight line, though this need not be the case). The figure 
indicates the projections 7 Q , 5 a and <5 r (A) of the cycles A a , B a and B r onto the z plane. 
The cycles A a can be identified with their projections j a . The cycles B a are defined such 
that, when crossing the cut I a going upwards along these cycles, one moves from the branch 
uq to the branch u± = —u$. Thus the 'lower halves' of these cycles lie on the branch uq, 
while their 'upper halves' lie on the branch u±. A similar convention is used for B r (A). With 
the orientation of the Riemann surface induced by its complex structure, this implies the 
intersections A a nB a = +1 (note that the cycles A a and B a intersect in a single point, which 
lies on the branch no). 

where B a are cycles on the large N Riemann surface chosen as explained in figure 6. 

Consider the cycles B a = J2 r p^a Bp for all a = 1 . . . r — 1. Then (4.32) implies: 

fxJP) = + I dzu(z) for a = 1 . . . r - 1 . (4.33) 



The quantity fi^ is undetermined and can be fixed arbitrarily. Following [1], we take 
jjif^ = § Br dzu(z), where A is a point close to infinity and B r (A) is the cycle described 
in figure 6. Defining B a (A) = B a + B r (A) for all a = 1 ... r — 1, we obtain: 

/^ 0) = U a for a = 1 . . . r , (4.34) 

with: 

U a := [ dzu(z) . (4.35) 



Relation (4.34) shows that the chemical potentials \x a are the finite A^ analogues of the 
periods ITq,. 

Note also that the filling fractions can be expressed as periods of the meromorphic 
differential udz over the cycles A a of figure 6: 

S a = i ^-.uJq(z) = I ^-.u (z) = I ^-.u(z) . (4.36) 

J^ a 2/K1 J^a LlXl JA a AlTl 

In the second equality, we used relation (3.22) and the fact that W(z) is a polynomial. 
Equation (4.15) now gives the special geometry relation of [1]: 

U a -U r = ^F . (4.37) 
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Note that the quantity in the right hand side is the planar limit of the microcanonical 
generating function. 

It is clear from figure 6 that A a fl B a = —A a n B a -i = +1 and A a fl Bp = if 
P ^ a, a — 1. This gives A a fl .Bg = — S a fl Ap = 5 a p. Since we also have A a (~) Ap = 
B a D Bp = 0, it follows that A a , Bp have canonical intersection form. Hence we have a 
canonical system of cycles A a , B a with a — 1 . . . r : 

A a n = -5 a n = 6 a>/3 , A a n Ap = B a n 5/3 = for all a = 1 . . . r . (4.38) 

This shows that the properties essential for the conjecture of [1] hold at the level of the 
holomorphic matrix model. 

5. Holomorphic ADE models 

In this section, we consider the case of holomorphic ADE models, focusing on the simple 
example of the holomorphic A 2 model. As mentioned in the introduction, the Hermitian 
approach to such models (and its attending regularization, discussed in [22]) leads to 
various technical problems which would violate the conjecture of [1]. The purpose of 
this section is to show explicitly how such issues are avoided in the holomorphic set-up, 
and to provide a reconstruction theorem similar to the one found for the one-matrix 
case. In particular, we shall extract explicitly the associated Riemann surface (which is 
expected from the work of of [3]) and show that one must use a certain 'renormalization' 
procedure in order eliminate unwanted constraints on its moduli. In fact, we shall find 
that the curve expected from the work of [3] and [23] can be obtained by using a 
certain regulator (which is natural in the holomorphic set-up), though only in the limit 
where this regulator is removed. As we shall see below, this limit can be described as 
a statistical ensemble of 'reduced' holomorphic models. 

5.1 Construction of the models 

Consider an ADE quiver diagram with nodes indexed by a = 1 . . . k where k is the 
rank of the associated simply-laced group. Consider N a x N a complex matrices $^ 
for each node, and N a x Np complex matrices Q^ a ^ for each pair of nodes a, (3 which 
are connected by an edge (in particular, we have two matrices Q( a ^ and for 
each edge of the quiver). We let s a p = sp a be the incidence matrix of the quiver and 
c a p = 25 a p — s a p be the associated Cartan matrix. 

By analogy with [22] , we define the holomorphic ADE matrix model associated with 
this quiver by: 




(5.1) 
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where W a are polynomials of degrees n a and: 

W int (*, Q) ■= E s ^ [tr(QW$WQW) - tr(Q^^Q^)} . (5.2) 

a</3 



The gauge group is: 



G:= ]jGL(N a ,C) , (5.3) 



with the the obvious action: 

Q W)) ^ (S a ^S-\S a Q^S^) (5.4) 

foiS a eGL{N a ,C). 

To completely specify the model, one must choose an appropriate integration man- 
ifold T. Before explaining our choice, let us comment on the Hermitian approach [22]. 
In that case, one takes T to consists of matrices $( a ), Q^ a ^ such that are Hermitian 
and = (Q^^Y for neighboring nodes a and (3. Such a prescription makes sense 

only if all W a have even degree (otherwise, the integral diverges because the absolute 
value of the integrand explodes when the norm of some is large). However, the 
Hermitian prescription immediately leads to other problems, arising from the integrals 
over Q. Indeed, it is easy to see that these will bring infinite contributions when some 
eigenvalue of coincides with some eigenvalue of for neighboring a 
and (3. In the Hermitian framework, the solution to this problem is to require that the 
eigenvalues of neighboring $( Q ) can never coincide — for example, by taking to 
have alternately negative and positive eigenvalues [22]. It turns out that this prescrip- 
tion would violate the conjecture of [1]. To understand why, consider for simplicity the 
Hermitian A 2 model (whose holomorphic version is studied below). The large N limit 
of this Hermitian model can be extracted by an argument which is formally identical 
to that presented in Appendix 3 5 , and is governed by an algebraic curve which is a 
triple cover of the complex plane. Denoting its three branches by u±, u 2 and u% (in a 
convenient enumeration), one finds that cuts connecting U\ and u 2 would correspond 
precisely to loci where equal eigenvalues of $W and $( 2 ) accumulate — a situation which 
is forbidden by the regularization prescription used to define the model ! Therefore, 
the regularization prescription of the Hermitian model requires that all such cuts are 
reduced to double points, which means that the large iV curve must always be singular. 
Moreover, this violates the requirements of the dual field theory, since it would require 
that some filling fractions vanish identically, thereby violating the study of gaugino 



5 Except that the third equation in (5.31) used in the appendix never plays any role for the regu- 
larization of [22]. 
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condensates performed in [23]. Of course, the nonzero cuts of the resulting curve would 
also be constrained to alternately lie on the positive and negative halves of the real 
axis. It should be clear from this discussion that Hermitian ADE models are quite 
unnatural for the conjecture of [1]. Below, we shall show how a certain regularization 
prescription of the holomorphic A 2 model allows one to obtain an ensemble whose large 
N limit satisfies the basic requirements of the conjecture. Unlike the Hermitian regu- 
larization of [22], the 'holomorphic' regulator used below can be removed in a manner 
which allows us to recover a smooth Riemann surface. 

Returning to the holomorphic model, we shall choose the integration manifold T 
as follows: 

(1) Fix contours 7 (a * ) in the complex plane such that: 

,y(<*) p| ryiP) — f or neighboring a and j3 . (5.5) 

and such that each 7^ connects two convergence sectors of W a (as defined in Section 
2). 

(2) Let A( 7 ) be the set of all matrices (D^ a \ QW>) (with a, (3 = \ ...k) which 
satisfy: 

(a) D^ = diag(AS a) . . . \<$), with distinct X^ . . . \%l 

(b) A^ lies along 7^ for each a and j 



3 



c )Q)i — i. a ) 3 (i3) (here the bar denotes complex conjugation). 

A ■ — A 



(3) Finally, we let V be the union of G-orbits 6 of elements of A under the action 
(5.4). Gauge- fixing the action (5.1) gives the eigenvalue representation: 

Z=[ <).../ dX^...[ d\[»\..[ d\%> e~ N ^L^ w ^ 

where we dropped constant prefactors and used the identity: 



duAdve^ x - x > v = -^- , (5.7) 

A — A 



with A 7^ A' complex and a the contour in C given by: 

« = xV (5 - 8) 

Condition (5.5) acts as a 'complex' regulator for the holomorphic model, by pre- 
venting common eigenvalues of neighboring and As in the Hermitian case, 

6 To be more precise, one has to take into account the action of an appropriate group of permutations 
on the set A. This can be done as in Appendix 1 by working with a fundamental domain of A under 
this discrete action. 
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working with the regularized model would therefore not suffice to recover the entire 
moduli space of planar solutions required by the conjecture of [1]. To eliminate the 
constraints, the regularization condition (5.5) must be removed by taking the limit 
of coinciding 7^. This limit can be performed in such a way that the end result is 
a 'renormalized' model which can be described as a statistical ensemble of 'reduced' 
holomorphic models. We now show how this works for the case of holomorphic A 2 
models. 

5.2 Example: the holomorphic A 2 model 

For an A 2 quiver, one has p = 2 and four matrices & 2 \ and Q( 21 \ The 

partition function takes the form: 

y 7 (D y 7 (D Wi y 7 ( 2 ) y 7 <*> N2 nv^-^f) 

(5.9) 

for two disjoint curves 7^ and 7^. 

5.2.1 Classical vacua 

Extremizing the action: 

S = NtriW^W) + W 2 (^)} + tr[Q<- 2 W>QM - qMqWqW] (5.10) 

gives the equations: 

NWl($W) = -Q (12) Q (21) , NW^) = QWQW , 
$( 2 )q( 21 ) = q( 21 )$« , $Wq( 12 ) = q( 12 )$( 2 ) . (5.11) 

Combining these, one easily obtains: 

w^'($ (1) )K($ (1) ) + w^'($ (1) )] = 

W^ 2 '($ (2) )K($ (2) ) + W^'($ (2) )] = . (5.12) 
Assuming that and $( 2 ) are diagonalizable with eigenvalues A^ and Xf\ one finds: 

^ 1 (aS 1) )K(aS 1) ) + ^(aS 1) )] = o 

W!>{Xf )[W[(Xf ) + W^Xf)] = . (5.13) 



On the other hand, the last row of equations in (5.11) gives: 

(Af - X?)Q^ = 

(A, (1) - Af = . (5.14) 
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In the generic case when W[ and W 2 have no common roots, equations (5.13) and 
(5.14) show that a typical classical vacuum is specified (up to permutations of indices) 
by choosing: 

(1) roots A?>...Aig_ fc oiW{, 

(2) roots \?...\$_ k of Wi, 

(3) roots Hi . . . Hk of W 7 ^ + W 2 such that A^_ fe+m = A^_ fc+m = At m for m = 1 . . . k. 

(4) some nonzero values for Q^-k+m^-k+m and Q&^. fc+roiJV2 _ fc+m for to = 1 ... A;. 
In fact, one can set Q^_ k+m)N , 2 _ k+m — 1 for all to = 1 ... k by using the gauge trans- 
formations. 

Less generic vacua arise, for example, by allowing some of the n k to coincide. 
Notice that vacua with k ^ (i.e. vacua for which some A^ coincide with some A^) 
are removed when imposing the condition 7^ fl 7^ = — such vacua will not be 
visible unless one removes this condition on the contours. We now proceed to remove 
this regulator, by taking the limit in which the two curves coincide. 

5.2.2 The limit of coinciding contours 

To study the limit when 7*- 1 - 1 and 7*- 2 - 1 coincide, we let 7 ( - 2 - ) = 7, 7*- 1 -* = 7 + r/n (where 
n is the normal vector field to 7, chosen as in figure 3 of Section 3) and take the 
positive quantity i] to zero (note that this requires 7 to asymptote to lines lying in 
the intersection of some convergence sectors of W\ and W<i). For small rj, we can use 
the length coordinate s of 7 = 7^ also as a parameter along 7^. We then have 
AW(s) = A(s) + rro(s), where A(s) is the parameterization of 7. For rj — > + , the 
Sokhotsky formulae (3.5) give: 

Aprhf = % s <VA( s f) " W) Ha ' ] ~ sf) ■ (5 ' 15) 

Therefore, the denominator in the integrand of (5.9) takes the following form: 



min(N lt N 2 ) r, (1) _ (2) x r, (1) _ (2) v 

E E A ff - k J -' n 

k = 1 < i! < ... < i k < AT X O-GEfc J • ' • A V A Jfc J (»,j)#(»l,j<r(l))-(U.Ja(fc)) 

1 < 31 < ■ ■ ■ < 3k < N 2 

+ ( higher incidence terms) , 

where is the group of permutations on k elements. The 'higher incidence terms' 
are terms involving delta-function products of the type 8(s^ — s^)5(s^ — s^) = 
5(s^ — s^ 2 ' ) )5(s^ 2 ' ) — sj^) with distinct i, j, k. Such terms can be neglected since — as 
we shall see in a moment — they do not contribute to the final result. 
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Substituting (5.16) in (5.9) gives: 

min(JVi,JV 2 ) AM /V ' 

where: 

Z fc ( 7 ) = I dm . . . I dfi k I dX? ... [ dX^_ k ( dXf ] ... [ dX%\_ k A k e~ NS * , (5.18) 
with: 

^ := A,P- 1 (1) (2) , (5.19) 

11 i= 1...JV! - k {\ ~ A j ) 
j = 1 . . . N 2 - k 



II fe-A? 1 ) II (Af'-ft) (5.20) 



i = 1 . . . JVi — fe i = 1 . . . N 2 - k 

j = 1 . . . k j — 1 ... k 



and: 

JVi-fe N 2 -k k 

S k := E W^XV) + E ) + E [^i(^) + ^)] • (5.21) 

j=i j=i j=i 

Notice that the 'higher incidence terms' of equation (5.16) bring zero contribution 
to (5.17). This is due to the extra-factors of /ij — fi k contributed by the two products 
in the numerator of the weight factor of (5.9). 

Observation Consider a quantity H(z) (which depends on A^ but is symmetric 
under separate permutations of A^ . . . A^j and of A^ . . . A^). Then its average (H(z)) 
has the following behavior in the limit 77 — > + : 



2^k=0 I / fe!(AT 1 -fe)!(iV 2 -fe)!^\- H \ Z )I 



k 



hm (H(z)) = . (M M , ' := (H(z)) lim , (5.22) 

2^k=o \ L1X ) fe!(Ari-fe)!(Ar 2 -fc)!^ fe w; 

where: 

(H(z)) k = f d/ii ... fdfM k f d\ { V ... f d\%l_ k f d\™ ... f d\%l_ k A k H(z)e- NS * . (5.23) 
Thus (H(z))ii m is simply a weighted average taken over the limiting ensemble. 
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5.2.3 Equations of motion for the limiting ensemble 

Let us fix a component k of the limiting ensemble, and work with the model defined 
by the partition function Z k . Writing as an exponential gives: 

Z k {i) =Jda 1 ...Jda k J ds? ...J d4l_ k J dsf ] ...J ds%l_ h e~ N ^ , (5.24) 

where: 

i = 1 . . . JVi - k i = 1 . . . N 2 - k i = 1 . . . iVi - k 

j — 1 . . . k j — 1 . . . k j — 1 . . . N 2 — k 

Ni — k N2 —k k 

-jfT, lnA ( s ? } ) -jrT, lnA (4 2) ) - ^E lnA ^) ( 5 - 25 ) 

j=i j=i j=i 

and we wrote A^ = A(s^),A^ = A(s^ 2 - ) ) and /ij = A(a"j). Extremizing this with 
respect to sf"^ and a { gives the equations of motion: 

2 JVl_fe ' 1 1 N2 ~ k 1 1 fc 1 



N 



y-v ' 1 1 ^ 1 1 

A^-AOf) " * £ A (s f))-A(4 2 )) + N £ A(fl (D) _ A(<7j) " 

= ^(A( S «))-1 X ^- 



k 1 1 iVi-fe , JV 2 -fc 



A(«, (1) ) 


-AOf) 






A( S f } ) 


-a(4 1} ) 


fc 


1 



^A( S «)2 



£ A( S f >) - A( S f) N £ A(s f)) _ A(a (D) N £ A(s f)) _ A(<Tj) ~ 

W »< A <'< )} ~iV A(s ( 2 ) )2 



- ^(A(a i )) + W^(A(a i ))- ' ' M ' T ' ! 



X(a t ) 

Let us introduce the spectral densities: 

P ia \s)--=^ k 8(8-8P) + ±j:5(8-* s ) , (5.27) 
with the normalization: 

r N 
J dsp^(s) = . (5.28) 

We also introduce the traced resolvents: 

1 Na ~ k 1 1 k 1 /■ r s"i 

^ W: =^ § 7^ + *£j=m ^ (5 - 29) 
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5.2.4 The large N Riemann surface 

For every quantity 0, we define the large N expansion coefficients 4>j as in equations 
(3.17). In particular, we have the planar limits and: 

4 a \z) = f ds^^ . (5.30) 
J z — \(s) 

In the large N limit, the eigenvalues A(sj^), A(s^) and X(cr k ) accumulate on curve 
segments sitting along 7. The planar limits (5.30) will have cuts along three types of 
curvilinear intervals: 

(1) loci Cf 3 resulting from a planar distribution of the eigenvalues X(sf^) with 
j — 1 . . . Ni — k; these will be cuts of ou^ 

(2) loci C 23 supporting a distribution of the eigenvalues A(s^) with j = 1 . . . N 2 — k; 
these give cuts of uo^ 

(3) loci resulting from a distribution of \(<Jj) with j — 1 . . . k] they give common 
cuts of uJq~^ and uof^ . 

Here a, b, c are indices counting the various occurrences of each type of cut. Note 
that the third type of locus can only arise from a component of the limiting ensemble 
for which has a finite limit as N — > 00. Thus cuts of type C 12 can only develop in 
the large N limit of the 'renormalized' model with = 7^. 

The planar limit of the equations of motion (5.26) gives: 

2 J ds '>w^h) -I ds 'xfih) = mMs)) ■ A(s) e CL (5 - 31) 



12 



These relations act as large iV saddle point equations for the limiting ensemble (5.17). 
They also represent the 'quantum' version of the three branches W{(\) = 0, W^A) = 
and W[(\) + W 2 {X) = of the classical equations of motion (5.13). It is also clear that 
we have: 

P ( o\s)=p£\s) , \(s)eC c 12 . (5.32) 

This is the analogue of the classical relations A^ = A^ on the corresponding branch 
of the moduli space (see point (3) of Subsection 5.1.1). 



31 



As explained in Appendix 3, one can use (5.31) and certain partial fraction decom- 
positions to derive the algebraic constraints: 



4\zf - u}? } {zH"(z) + ^\zf - W[{z)^\z) - W>{z)^\z) + ft\z) + tf'{z) = 



,( 2 ) 



,(2), 



r(2), 



4 1 \z) 2 4 2 \z)-Wi(z) ^\zf + f^\z)-W[{z)^\z) +^(«)-(1^2) = (5.33) 



where /o Q \<?o^ are polynomials defined through: 



(«)/ 



dsp{, a) (s) 



,^)-^(A(«)) 



z-X(s) 



and 



^(z) := | da/ dp[W{(z) - W[{\{a))} 
g!?\z) := fdaf d(3\Wfc) - W^\(a))] 



pt\<x)pi\p) 



(A(a) ->(/?))(* -A(a)) 
Note that the normalization conditions (5.28) imply: 

lcoeff(/ (a) ) = -^deg(WQlcoeff (W a ) . 
where lcoeff(...) denotes the leading coefficient. Writing: 



where: 



u$-\z) 



h{z) :-- 
h{z) :-- 



= ui{z) + ti{z) 
= -u 2 (z) + t 2 (z) 

2Wj(z) + W' 2 (z) 
3 

2W 2 '(z) + W[{z) 



brings the constraints (5.33) to the form: 

ui(z) 2 + ui(z)u 2 {z) + u 2 (z) 2 = p(z) 
ui(z) 2 u 2 (z) +u 1 (z)u 2 (z) 2 — —q(z) 



where: 



(5.34) 



(5.35) 



(5.36) 



(5.37) 



(5.38) 



p(z) := tl (z) 2 - h{z)t 2 {z) + t 2 (z) 2 - ff \z) - f i2 \z) 
q{z) := -h{z)t 2 {z) [t 1 {z)-t 2 {z)]-t 1 {z)f^{z)+t 2 {z)fi l \z)-g,{z) . 



(5.39) 



(5.40) 



32 



In the last equation, we introduced the polynomial go(z) = g^\z) — g^\z) which has 
the following explicit form in terms of matrix model data: 

g (z) = [da [d(3[W[(z)-Wi(X(a))\ " ! ' (n)A<r ' ( ° 



(X(a) - X(P))(z - X(a)) 
-J da J dp W(z) - WM*))\ (A(a) ^ ( gC ( -A(.)) (5 ' 41) 

Defining 1x3(2;) := —1x1(2) — u 2 (z), identities (5.39) can be recognized as the Viete 
relations of the cubic: 

u 3 - p{z)u - q{z) = (5.42) 

when the left hand side is viewed as a polynomial in u. This shows that 1x1(2), 1x2(2) 
and 1x3(2) are the three branches of the affme algebraic curve (5.42). This is the precise 
form of the curve suggested in [3] (where the explicit form of the polynomials p, q in 
terms of matrix model data was not given). Note that the left hand side of (5.42) can 
also be written as: 

(u- Ul )(u-u 2 )(u-u 3 ) = (u+tJiu-hjiu-h+^ + ifW + / (2 V+*i/o (2) -* 2 /o (1) + <?0 • 

(5.43) 

Studying the behavior of Uj near the branch cuts of (5.42) allows one to identify 
these as the loci Ci 3 , C 23 and C\ 2 where the eigenvalues accumulate; then the jump 
equations across these cuts can be seen to be equivalent with the planar equations of 
motion (5.31). Below, we shall use this reasoning in order to give a reconstruction 
theorem for the holomorphic A 2 model, similar to the one we found in Section 3 for 
the holomorphic one-matrix model. 

5.2.5 The reconstruction theorem 

Let us start with a curve of form (5.42) with p, q given by (5.40), where and g 
are complex polynomials of degree n — 2 subject to the constraints (5.36). Given such 
data, we show how one can construct a curve 7 and distributions Pq"\s) along this 
curve such that relations (5.28), (5.30), (5.31), (5.32) and (5.34), (5.41) hold. 

Using expression (5.43), one finds the following asymptotic behavior for large z: 

de^Wi)lcoeff(Vyi) 

u,{z) = +t 2 + - (Z e ^l] w , + 0(1/ z^) . (5.44) 
deg(W 2 )\coeS(W 2 ) 
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In particular, we can use these asymptotic forms in order to fix the indexing of branches 
for (3.23) (i.e. U\ is the branch which asymptotes to —t±, while u 2 is the branch which 
asymptotes to +t 2 ). 

When viewed as a triple cover of the z-plane, the curve (5.42) has three types of 
cuts, namely those connecting the pairs of branches (1,3), (2,3) and (1,2). Denote 



these cuts by C" 3 , C b 23 and C{ 2 . Picking 7 to contain all cuts, we let 00^ and u^' be 
given by relations (5.37) and define: 

, for X(s) e C? 3 
for X(s) e C b 23 , 

(5.45) 
for X(s) G Cfa . 



,( 2 ) 









-iO) 


-4 1} (A(s)+^0) 






;4 2) (A( S ) 


-iO) 


-4 2) (A(s)+iO) 




■= — A(s) 


;4 X) (A( S ) 


-iO) 


-4 1) (A(s)+iO) 




= — A(s) 
2ni V ; 


;4 2) (A( S ) 


-iO) 


-4 2) (A(s)+iO) 



The very last equality in these relations follows by using (5.37) and the fact that 
u 3 = —Ui—u 2 is continuous across cuts of type Ci 2 . In particular, this means that (5.32) 
holds. Extending p^f 1 by zero to the entire curve 7, relations (5.45) and the Sokhotsky 
formulae immediately imply that equations (5.30) hold. In turn, these relations, the 
asymptotic behavior (5.44) and the constraints (5.36) on the leading coefficients imply 
the normalization conditions (5.28). 
We next notice that: 



ui{z) - u 3 {z) 
u 2 {z) - u 3 (z) 
ui(z) - u 2 {z) 



(1), 



2u 1 {z) +u 2 {z) = 2u y o'{z) - Jg\z) - W[{z) 
Ul {z) + 2u 2 {z) = u£\z) - 2uf\z) + W&z) 
4 1 \z)+4 2 \z)-(W{(z) + W 2 '(z)) , 



(5.46) 



where we used relations (5.37) and (5.38). Since a cut of type C a p connects the branches 
u a and up, we have the jump equations: 

u a (\ ± 0n(A)) = up(X =F iOn(X)) =>• (5.47) 
[u a (X + 0n(A)) - up(X + 0n(A))] + [u a (X - 0n(A)) - up(X - 0n(A))] = 

along such a cut. Combining this with relations (5.46) and using equations (5.30) and 
the Sokhotsky formulae (3.5) shows that satisfy the planar equations of motion 
(5.31). 

To prove that (5.34) and (5.41) hold, one can now simply repeat the original argu- 
ment (found in Appendix 3) leading to the large N curve (5.42); this is possible since we 
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already showed that all assumptions of that argument (namely relations (5.30), (5.31), 
(5.32)) hold. This shows that equations (5.39) and (5.40) must hold with polynomials 
/, g defined by relations (5.34) and (5.41). Since by hypothesis the same relations hold 
with the original polynomials flf 1 and g , it immediately follows that the latter can 
indeed be expressed in the form (5.34) and (5.41). This concludes the proof that the 
'renormalized' holomorphic A 2 model indeed probes the whole moduli space of (5.42) 
with the constraints (5.36). 

Observation The third 'branch' of the planar equations of motion (5.31) is only 
allowed in the 'renormalized' model described by the limiting ensemble (5.17). If one 
works with the regularized model instead (the original model for which 7^ and 7^ 
are disjoint), then cuts of type C{ 2 are not allowed. Indeed, such cuts connect the 
branches U\ and u 2 , and thus are cuts for both uj^ and uj^\ which would require 
C\ 2 C 7*- 1 - 1 H in contradiction with the regularization condition 7^ n 7*- 2 - 1 = 0. 
Hence the regularized model can only probe that part of the moduli space of (5.42) for 
which all cuts of type C\ 2 are reduced to double points. This is similar to what happens 
for the Hermitian A 2 model, as discussed at the beginning of this section. In particular, 
the regularized representation (5.9) cannot capture the entire family of curves (5.42), 
and, in fact, can only describe a subvariety in the space of all allowed polynomials 
and go*\ To explore the full moduli space, one must consider the limiting ensemble 
(5.17). In the context of the Dijkgraaf-Vafa conjecture, the dual field theory is an 
N — 1 supersymmetric field theory derived from an A 2 quiver (such theories have been 
studied in [23, 24]). This field theory is certainly allowed to explore the branch whose 
classical description is given by W[(\) + W 2 (X) = 0. Therefore, consideration of the 
limiting ensemble (5.17) of the holomorphic model is required by the conjecture of [1] 
as applied to the A 2 quiver field theories of [23]. A similar analysis can be performed 
for general ADE models, with the same conclusion. 
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A. Integration over gauge group orbits for the holomorphic 
one-matrix model 

In this Appendix, we show how the eigenvalue representation (1.14) results from the 
matrix integral (1.7). 
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A.l Orbit decomposition of M. 

Let A be the set of diagonal N x N complex matrices D with distinct eigenvalues: 

A := {D = diag(Ai . . . \ N )\>^j G C and J] (A; - Aj) + 0} (A.l) 

and fix a fundamental domain A for the obvious action of the permutation group Sjv 
on A. Then relation (1.2) gives the orbit decomposition: 

M = U DeAo O D , (A.2) 

where Od is the GL(N, C)-orbit of a matrix D G A under the similarity action (1.8). 

The stabilizer of each D G A in GL(N, C) is the subgroup T N ps (C*)^ consisting 
of diagonal matrices. Hence each orbit has the form: 

O d = H:=GL(N,C)/T n , (A.3) 

where the homogeneous space H has dimension A^ 2 — N (here Tjy acts on GL(N, C) 
from the right, i.e. S — > ST for S G GL(N, C) and T G T N ). The orbit decomposition 
(A.2) takes the form: 

M = H x A . (A.4) 

A.2 Decomposition of w 

Let: 

N 

w A = Y[ d\j (A. 5) 

be the translation-invariant holomorphic volume form on A and: 

w H = K^jUij (A. 6) 

be the left- invariant holomorphic volume form on the complex homogeneous space H, 
where uj = S^dS is the matrix whose elements give a basis of left-invariant holomorphic 
one-forms on GL(N, C). 

If ls{S') '■= SS' is the left action of GL(N, C) on H, then wh satisfies: 

(Is)*wh = w h , SeGL(N,C) . (A.7) 

Using the projections n H and n A of M. = H x A onto its two factors, we define a 
holomorphic A^ 2 -form on M. by: 

w := ir* H (w H ) A 7r^(w A ) . (A. 8) 
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For dimension reasons, this must be related to w through: 

w = fw , (A.9) 

where f{M) is a holomorphic function on M.. Using the left-invariance (A. 7) of wh, it 
is easy to check that wo is invariant 7 under the action (1.8): 

r(S)*(w )=w . (A.10) 

Using GL(N, C) invariance of w and w , relation (A.9) implies that / is invariant: 

f(SMS- 1 ) = f(M) . (A.ll) 

In particular, we have f(M) = f(D) if SMS^ 1 = D with D diagonal. It thus suffices 
to determine the values of / for diagonal matrices D. 

For this, we first describe the the cotangent space to M. at the points of A (where 
A is viewed as a subset of M). In the vicinity of A , we can write M = SDS^ 1 pa 
{l + A)D{\ — A) pa D + [A, D], where S = e A pa 1 + A with A an infinitesimal generator 
of GL(N, C). Therefore, at a point D of A , we have: 

dM = dD + [dA, D] dM^ = 5 i3 d\i + (Xj - X^dA^ . (A.12) 

Note that there is no [A, dD] piece in the right hand side of this relation, since we 
compute the form dM at the point D (where A = 0). 

Noticing that wh = A^jdAij at such a point, relation (A.12) gives the form of w 
at M = D : 



w 



= (_l ) ^-D/2 Aij . dMy = ( _ 1) 



7V 2 (iV-l)/2 



n - x ^ 



w . (A.13) 



Comparing with (A.9), we find that / is given by the usual Vandermonde determinant: 

/(M) = /(£>) = J] - A,-) . (A.14) 

Combining this with (A.9) gives the following expression for the holomorphic measure: 

w = (-l) N ^ N - 1 ^H(X t -X J )w , (A.15) 

where Aj are the eigenvalues of the matrix M at which we evaluate w. 



7 To check this, notice that tth ° t(S) = Is ° t^h- 
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A. 3 The eigenvalue representation 

The GL(N,C) invariance of the action, relation (A. 15) and the decomposition (A. 8) 
allow us to perform the integral over the gauge-group variables in the partition function 
(1.7): 

1 r N 

z N ( % t) = u-ir^-^hvoi{ H ) n ^ n - ^ n ^ wm (a.i 6 ) 

where hvol(H) = J h u>h is the 'holomorphic volume' of if and: 

A ( 7 ) = {D = diag(A! . . . X N ) e A |A j e 7 , V? = 1 . . . N} . (A.17) 

Since this relation holds for any choice of fundamental domain A , we can write (A. 16) 
in the form: 

Z N ln,t) = ^{-l^W-W-L hvol{H)Z N {^t) , (A.18) 

where: 

r N x^n 

zs{%t)= / n^n^-A^^-^ (A.i9) 

J Hi) j= i #j 

with: 

A( 7 ) = {D = diag(Ax . . . X N )\\j G 7 , Vj = 1 . . .N and J| (A ; - \) + 0} (A.20) 

This gives the eigenvalue representation (1.14). 

Observation When writing the representation (1.14), we have tacitly extended the 
integral to allow for colliding eigenvalues \ = Xj (this is certainly allowed, since the 
integrand of (A.19) is well-behaved there). This amounts to treating non-diagonalizable 
matrices by a point-splitting prescription, and can be formulated in more detail by 
working with: 

A e := {D = diag(Ai . . . X N )\\j G C and |A; - Aj| > e for i ^ j} (A.21) 

instead of A and with a similar modification }A e of M.. Then one defines the partition 
function as the limit e — > + of the regularized partition function obtained in this 
manner. It is easy to adapt the derivation given above in order to include explicitly 
such a regulator. The result, however, is the same as (1.14), because the integrand of 
(A.19) is well-behaved when eigenvalues come close to each other. 
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B. Example of the relevance of convergence sectors: the case 
of a cubic potential 

Consider the potential: 

W(z) =t 3 z 3 + t 2 z 2 . (B.l) 

This example appeared in the paper [6], where it was used to carry out a one-loop test 
of the Dijkgraaf-Vafa conjecture. As in [6], we assume for simplicity that t 2 and t 3 are 
real and positive, and write them as t 2 = y and £3=3 with positive m and g. The 
potential has two local extrema along the real axis, namely a local minimum at a\ — 
and a local maximum at a 2 = — ||| = — — < 0. Also note that W(z) approaches ±00 
as z approaches ±00 along the real axis. 

B.l Summary of the procedure of [6] 

The paper [6] follows [1] by formulating a conjecture mapping a one-matrix model based 
on the potential (B.l) to the closed topological B-model on a certain non-compact 
Calabi-Yau space. When testing this, [6] performs a perturbative expansion of this 
model around the two extrema a\ and a 2 , thereby re- writing the theory as a a two- 
matrix model for matrices M 1 ,M 2 with cubic potentials Wi(Mi) and W 2 (M 2 ) and a 
logarithmic interaction term Wi n t(Mi, M 2 ). Since the point a 2 is a local maximum 
for the original potential W, this would produce a negative definite quadratic term in 
W 2 (M 2 ), if one starts with a Hermitian matrix model, thereby leading to an ill-defined 
perturbation expansion. To cure this problem, [6] proposes that one should take M x 
to be anti-Hermitian and M 2 to be Hermitian. With this redefinition, the authors 
of [6] compute the first few perturbative terms and find agreement with the one-loop 
computation of the Kodaira-Spencer theory of the Calabi-Yau dual. We now show how 
this pragmatic procedure can be recovered in the holomorphic framework. 

B.2 Justification in terms of holomorphic matrix models 

Remember that a Hermitian matrix model based on the cubic potential (B.l) is ill- 
defined. This is due to the exponential increase of the integrand for Xj — > —00. On 
the other hand, the holomorphic matrix model leads to a meaningful integral, provided 
that one chooses appropriate asymptotic sectors. 

Note that our definition (1.7) uses the weight Y\i^j(Xi — Xj) rather than the weight 
Tii<j(Xi — Xj) 2 which is used in [6]. In this appendix, we shall temporarily switch to 
the representation: 

Z N (>y,t) = I dX l ... I dX N H(X t - Xjfe- N ^U w ^ (B.2) 

J ■y J 7 - , - 
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which differs from our convention (1.14) by a sign prefactor of (— l) Ar ( Ar ~ 1 )/ 2 . This 
avoids some annoying sign factors when comparing with [6]. 

Thus we start with the partition function (B.2) for the potential (B.l). Since 
£3 > 0, this model has #3 = and the asymptotic sectors shown in figure 7. We claim 
that the correct partition function relevant for the work of [6] is Z(l,0,t), associated 
with the asymptotic sectors k- = 1 and k + = 0. This is convergent by the analysis of 
Section 2. 

In this set-up, the prescription of [6] can be recovered as follows (we shall neglect the 
gauge group volume prefactors, which are irrelevant here). Since the partition function 
(B.2) depends only on the asymptotic sectors of 7 (namely z/_ 6 A\ and u + G A ), 
we are free to choose this curve as shown in figure 7. This asymptotes to some line 
d- with tangent z/_ e A\ for t — > —00, then follows a vertical line through the point 
a 2 = —m/g, and finally loops back to touch (and then follow) the real axis at some 
point x lying in between a 2 = —ni/g and a x = 0. Thus an eigenvalue A sitting on 7 
will be imaginary if it lies close to a 2 and real if it lies close to ao, thereby naturally 
implementing the requirement of [6]. The curve segments along 7 around the points a\ 
and a 2 for which A has these properties can be maximized by making the 'well' at the 
bottom of this curve to be very thin (i.e. take x to be very close to a 2 ) and very deep. 





t 
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Ql = » 

B 2 
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Figure 7: Convergence sectors for the case deg W = 3, #3 = and a good choice of contour 
belonging to the sector k + ) = (1,0). 

Following [6], we now expand the integral (B.2) (with 7 chosen as above) around 
the saddle point configurations: 

Aj 0) = a x = for i = l...N 1 , \f = a 2 = -m/g for j = N x + l...N (B.3) 
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with Ni some integer in the set {1 . . . iV}. Let N 2 := N — N\. This gives: 



Z & 7) = E TTTTTT I d^.-.f dfi Nl fdu^.J du Na A(fi, u)e~ NS (B.4) 

N 1 +N 2 =N iV l !iV 2l ^7 J 7 ^7 " , 7 

where we wrote Aj = A^ + fa, X^+j = ^iv|+j + u j an d ; 

iVi iv 2 m 

A (^ z/ )= n (^i-^ 2 ) 2 n k - ^- 2 ) 2 n n - ^ + — ) 2 

l<ii<i 2 <A r i l<ji<J2<iVi i=lj=l 9 

Ni N 2 3 

5 = E^)-E w (-^) + ^2 • (B.5) 

i=i j=i °y 

Treating /xj and z/j as small fluctuations, we naturally have /Xj G ?R and z/j G R, 
since the eigenvalues A lying along 7 are imaginary in the vicinity of a 2 and real in the 
vicinity of a\. Writing the last term in A as an exponential, one obtains a logarithmic 
interaction potential 8 : 



W int = 2N,N 2 In - + 2 £ £ In 1 + ^ - i/,- 



5 



i=l i=l 



(B.6) 



and one-matrix potentials W^i(/i) = W(fx) and W 2 (u) = — W(— v), which allows one to 
write the result as an interacting two-matrix model [6]: 

Z = V m ( dM x f dM2e -Ntr[W 1 (M 1 )+W2(M2)]+W mt (M 1 ,M2) 

N 1+ N 2 =N Ni\N 2 \ JM+=-M 1 Jm+=+M 2 

(B.7) 

where the first integral is over anti-Hermitian matrices while the second integral is 
over Hermitian matrices. Note that anti-Hermiticity of Mi arises naturally in the 
holomorphic set-up. This gives a conceptual justification for the procedure of [6]. Note 
also that we have provided a non-perturbative construction of the model involved in 
that work: it is simply the holomorphic matrix model with potential (B.l), considered 
in the 'phase' (k-,k+) = (1,0). 



C. Derivation of the planar constraints for the A2 model 

Let us show how the non-hyperelliptic Riemann surface expected from the observations 
of [3] arises in the holomorphic A 2 model. For this, we derive two algebraic constraints 

8 The authors of [6] further expand the logarithm as a power series in ^(/Xj — Vj), a procedure which 
is justified if - Vj)\ < 1. 
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which hold in the planar limit, as a consequence of the planar equations of motion 
(5.31). 

To derive the equations of interest, we shall use the partial fraction decompositions: 



(z — u)(z — v) 



u — V 



z — u z — V 



(C.l) 



and: 



(z — u) (z — v) (z — w) (u — v)(u — w) z — u (v — u) (v — w) z — v (w — u)(w — v) z — w 

(C.2) 

C.l The first constraint 

Using (C.l), one finds: 



^\zf 


= 2 J ds J ds 




= 2 j ds J ds 




= J ds J ds' 



\{s) - A(s') z-X(s) 
p { o\s)p ( i\s') 1 
X(s)-X(s') z-X(s) 



(C.3) 



A* 1; (*KV) 



+ 



pS^UV) 1 



X(s) - X(s') z - X(s) X(s) - X(s') z - X(s) 
Combining these equations gives: 

4%f - + W 2 = / ^"M^f + / *ft w «7=^f 

(C.4) 

To arrive at this relation, we decomposed the integrals over ds in (C.3) into the pieces 
corresponding to the cuts C" 3 , and C\ 2 . Then equation (C.4) results upon com- 
bining these pieces appropriately and performing the s' integral by using the planar 
equations of motion (5.31) and relations (5.32). We next write (C.4) in the form: 

J£\z? -J£\z)J£\z)+^ = , 

(C.5) 

where we used the planar equations of motion (5.31) and the definition of ivq*\z) and 
we introduced the polynomials: 



fi 2 \z) := / dsp^(s) 



.(2), 



X(s) 

Wj( z)-W&\(*)) 
z-X(s) 



(C6) 
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C.2 The second constraint 

To derive the second constraint, we use (C.2) to compute: 

(i), ^2 (2), x f , I AR f , 2p<, 1) (a)p( 1) (/j)p( 2) ( 7 ) 

4 >(*) 4 V) = y ^y ^y ^ (A(a) _ A(/3))(A(a) _ A(7))(z 



(A(a)-A(0))(A(a)-A( 7 ))(*-A(a)) 
+ jdajd(3j *7 , x(a) _ M0))(Ma) _ x , l))(z _ x(a)) > (C7) 



(A(a)-A(/?))(A(a)-A( 7 ))(*-A(a)) 

where we used redefinitions of 7) by permutations to bring the right hand side 
to a convenient form. Combining this with the equation obtained by permuting the 
indices 1^2 gives: 

4WW-#>(.H W M' - / rfa /^(A (a) ) (A( /Wg^ (a)) -(i ~ 2) , 

(C.8) 

where we used the planar equations of motion (5.31) to perform the integral over (3. 
Defining the polynomials: 

IX) („\ ._ / a~ / Aamr'(,\ _ W(\(r,\w Po \ a )Po \P) 



9b \ z 



J da J dp [W[{z] - W[{\{a))\ 



(A(a)-AG9))(z-A(a)) ' 
f%\<*)$(P) 



^(z) := fdafwWiz) - ^(A(a))] ^ ^ ( ^ A(a)) , (C9) 



we find: 



4\zfujf\z) - 4\z)^\zf + $\z) - g<?\z) - W'^U^z) + W>(z)U 2 (z) = , 

(CIO) 

where: 

Po 1) («)Po 2) (/?) 



Ul(z) '= J da J d(3 (X(a) - m)(z - A(a)) 

Using the equations of motion (5.31), these quantities can be written: 

U x {z) = f£\z) - Wi(z)4 1 \z)+4 1 \z) 2 

U 2 (z) = tf\z) - Wi(z)4 2) (z) + uf\zf . (C.12) 
Therefore, equation (C.10) becomes: 

4 1) (z) 2 4 2) (^)-w 1 (z)[4 1) (^) 2 -w 1 / (^)4 1) (^)+/o (1) (^)]+^ 1) (^)-(i - 2) = . (c.i3) 
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